Accuracy of direct gradient sensing by cell-surface receptors 
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' Chemotactic cells of eukaryotic organisms are able to accurately sense shallow 

chemical concentration gradients using cell-surface receptors. This sensing ability 
is remarkable as cells must be able to spatially resolve small fractional differences 
in the numbers of particles randomly arriving at cell-surface receptors by diffusion. 
An additional challenge and source of uncertainty is that particles, once bound and 
released, may rebind the same or a different receptor, which adds to noise without 
providing any new information about the environment. We recently derived the 

■ fundamental physical limits of gradient sensing using a simple spherical-cell model, 
y^Q . but not including explicit particle-receptor kinetics. Here, we use a method based 

on the fluctuation-dissipation theorem (FDT) to calculate the accuracy of gradient 
sensing by realistic receptors. We derive analytical results for two receptors, as well 

■ as two coaxial rings of receptors, e.g. one at each cell pole. For realistic receptors, 
we find that particle rebinding lowers the accuracy of gradient sensing, in line with 
our previous results. 

X 

I. INTRODUCTION 

Cells are able to sense gradients of chemical concentration with extremely high sensitivity 
and accuracy. This is done either directly, by measuring spatial gradients across the cell 
diameter, or indirectly, by temporally sensing gradients while moving [lj|. In temporal 
sensing, a cell modifies its swimming behavior according to whether a chemical concentration 
is rising or falling in time 0. This mode of sensing is typical of small, fast moving bacteria 
such as Escherichia coli, which can respond to changes in concentration as low as 3.2 nM 
of the attractant aspartate In contrast, direct spatial sensing is prevalent among larger, 
single-celled eukaryotic organisms such as the slime mold Dictyostelium discoideum (Dicty) 
and the budding yeast Saccharomyces cerevisiae j3, E|. Dicty cells are able to sense a 
concentration difference of only 1-5% across the cell |6(, corresponding to a difference in 
receptor occupancy between front and back of only 5 receptors [7| • Spatial sensing is also 
performed by cells of the immune system including neutrophils and lymphocytes jl] , as well 
as by growing synaptic cells and tumor cells. Interestingly, a direct spatial mode of sensing 
has also been demonstrated for the large oxygen-sensing bacterium Thiovulum majus 0j, 
indicating that direct gradient sensing is widespread among the different kingdoms of life. 



2 



There has been great progress in understan ding the limits of concentration sensing in 
bacteria such as E. coli, pioneered by Berg 0, [Tol and in understandin g th e origins 



of sensitivity in the underlying signaling network, pioneered by Bray [12|, [13|, [14|, [15|, and 



followed by others [16l . 1 171 . Il8l . Il9l . |20| . l2ll . |22J. By contrast, very little is known about what 
determines the accuracy of direct gradient sensing by eukaryotic cells. Recently, we derived 
the fundamental physical limits of direct gradient sensing, where the accuracy is limited by 
the random arrival of particles at the cell surface due to diffusion 23[. We used as models 
a perfectly absorbing sphere and a perfectly monitoring sphere (a la Berg and Purcell (llj). 
In these two models, gradients are inferred from the positions of particles absorbed on the 
surface of a sphere or the positions of freely diffusing particles inside a spherical volume, 
respectively. The latter case simulates rebinding of particles, as particles can enter and exit 
the spherical volume freely. In comparison, for the perfectly absorbing sphere, previously 
observed particles are never remeasured. As a result, we found that the perfectly absorbing 
sphere is superior to the perfectly monitoring sphere, both for concentration and gradient 
sensing (Table I). 

The superiority of the absorbing sphere may help explain the presence at the surfaces of 
cells of signal degrading enzymes, such as PDE for cAMP in Dicty and BAR1 for mating 
factor a in S. cerevisiae. Those surface enzymes could reduce or eliminate rebinding (and 
therefore remeasurement) of the same signal molecule. Quantitatively, our theory compares 
favorably to recent measurements of Dicty moving up shallow cAMP gradients 0], suggesting 
that these cells operate near the physical limits of gradient detection. 

While our recent models of the absorbing and monitoring spheres allowed us to derive 
the fundamental limit of gradient sensing, the models neglect the details of biochemical 
reactions, such as particle-receptor binding and downstream signaling events, which might 
further increase measurement uncertainty. To study the effects of particle-receptor binding, 
we here extend a formalism for the uncertainty of concentration sensing recently developed 



by Bialek and Setayeshgar [22], [24 



to the case of gradient sensing. This formalism uses 
the fluctuation-dissipation theorem to infer the fluctuations of the receptor occupancy (and 
hence the accuracy of concentration sensing) from the linear response of the average receptor 
occupancy to changes in receptor binding free energies. The effect of particle rebinding 
is included by coupling particle- receptor binding to the diffusion equation [24| . leading to 
correlations in time among the receptors. We report analytical results for two receptors (Fig. 
1), as well as two coaxial rings of receptors, e.g. one at each cell pole (Fig. 2). By assuming 
diffusion-limited particle binding to the receptors, we are able to directly compare to our 
previous model for the fundamental limits of gradient sensing. For realistic receptors, we 
find that particle rebinding lowers the accuracy of gradient sensing in line with our previous 
results for the absorbing and monitoring spheres (Table I). 



II. METHODS 



Bialek and Setayeshgar 24 1 previously presented a method based on the fluctuation- 
dissipation theorem (FDT) 25| to calculate the accuracy of measurement of chemical con- 
centration by receptors. We extend their method to calculate the accuracy of measurement 
of concentration gradients, and derive analytical results for (1) two receptors and (2) two 
coaxial rings of receptors, e.g. one at each cell pole. We start our derivation by considering 
an arbitrary number m of receptors. 

The kinetics of the ensemble-average occupancy rij(t) of receptor j due to binding and 
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Measurement uncertainty 



Perfect absorber 



Perfect monitor 



Ratio absorber/monitor 



Concentration: 



Gradient: ^fQ- 

{c /a) 2 



iirDac T 
1 

4irDac T 



5nDac T 
15 



7wDac T 



12. 

5 " 
60, 

7 ' 



2.4 
8.6 



TABLE I: Uncertainties in measured concentration and concentration gradient for two idealized 
cell models: a perfectly absorbing sphere (second column) and a perfectly monitoring sphere (third 
column). Also provided is the ratio of the uncertainties of the absorber and monitor. Parameters: 
diffusion constant D, radius of sphere a, averaging time T, and average chemical concentration cq. 
Table reproduced from Ref. 



23]. 




FIG. 1: Model for gradient sensing by two individual cell-surface receptors. Particles diffuse 
through the medium and randomly bind to and unbind from receptors, e.g. see the two sample 
paths indicated by black arrows. Particles may rebind the same (1) or a different (2) receptor. The 
receptor separation is given by r. 



unbinding of chemical ligands at the local concentration c(xj,t) is given by 

dnj(t) 



dt 



k + c(xj,t)[l — rij(t)] — k_rij(t). (1) 



Linearization about the mean steady-state occupancy rij = k + cj/ (k+Cj+kJ) at concentration 
Cj = c(xj , t) gives 

d{5rij(t)) = -(k + Cj + k_)5rij + cj(l - fij)8k + - nj5k_ + k+(l - n^Scj. (2) 

(Jjv 

By thermodynamics, the ratio of binding and unbinding rates is related to the free-energy 
difference Fj between the unbound and bound states of the receptor according to 

^-expfiV (3) 



k- yksT f 

Variations of the rate constants in Eq. [2] are equivalent to a variation or external perturbation 
of this receptor free-energy difference 

5F, 5k + 5k„ 

— J - = — . (4) 

ksT k + k_ 



Combining Eqs. [2] and HI one obtains after Fourier transforming in time to obtain a frequency 
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representation 



— iuo5hj(uo) = —(k + Cj + k-)6hj(uj) H — — - — =p — -5Fj(u) + k + (l — nj)5cj(u), (5) 

k B T 

where, e.g., 5hj{uj) describes the variations in the receptor occupancy as a function of fre- 
quency u). The dependence on 5cj(u) can be eliminated by using the (linearized) diffusion 
equation [24| in Eq. \5\ 



Of 

Using the Fourier transforms 

5c(x, t) 
5(x-xi) 



d(5c(x, t)) = DV , 5c ^ t) _£ 6{s _ £i) d{5 ni {t)) 



1=1 



dt 



(6) 



du 



d 3 k 



2tt J (2tt) 3 
d 3 k 



ik(x—xi)—k/A 



(2tt) 3 
Zn 



(7) 
(8) 
(9) 



where k = \k\ and where we have introduced a convergence factor A > in Eq. [8] to regulate 
the S function (effectively assigning a size scale ~ 1/A to the receptor), Eq. [6] yields 



Sc(cu, k) 



luJ 



Dk 2 — iuj 



i=i 



-ikii 



-k/Ag 



(10) 



Inverting the spatial Fourier transform back into real space, one obtains the concentration 
fluctuations at the locations of the receptors in terms of the occupancy fluctuations 



5c(xj,u>) = iu 5hi(u 
:=i 



^3^, e ik(xj-xi)-k/A 



(2vr) 3 Dk 2 -iu 



IU 



d 3 k 



-k/A 



(2tt) 3 Dk 2 - iu 



0J 



IU 

2^2 



k 2 dk 
Dk 2 — iu 



^3fc e ik(xj-xi) 

(2vr) 3 Dk 2 - iu 

k sm(k\xj — xi\) 
Dk 2 — iu 



dk 



(11) 
(12) 

(13) 



The cut-off A = tt/s, which accounts for the physical dimensions of the receptor, is used in 
Eq. [13] to set an upper limit of integration [24| . 

Following Bialek and Setayeshgar 24| , we imagine that the mechanism that reads out the 
receptor occupancy averages over a time r long compared to the correlation time between 
binding and unbinding events of a receptor (see below). In this case, we can apply the 
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FIG. 2: Model for gradient sensing by two polar rings of surface receptors. Particles diffuse through 
the medium and randomly bind to and unbind from receptors. Particles may rebind the same or 
different receptors (e.g., path of black arrows). Rings of receptors are separated by distance r. 
Also shown is receptor dimension s. 



low-frequency limit (u « DA 2 ). Using the expression 



lim 



±271-5(0), 



where v is real and P is the principal value of the associated integral, Eq. [13] becomes 



(14) 



5c{xj,u) = 7r^ 5n i( u > 



2n 2 D 

Inserting Eq. [15] for 5cj (x, u) in Eq. [5] yields 

iuoA 



E Shijco ) 



\Xj Xn 



(15) 



— iuo5hn 



k+(l — n 



jJ 2tt 2 D 



-/C_i_fl — Ua 



III! 



A-kD 



(k + Cj + k-) 
Shi 



5hn 



k B T 5F - 



(16) 



which depends only on the 5hj(u) and their conjugate variables 8Fj(u). These m equations 
describe how deterministic frequency-dependent changes in the free-energy differences 
5Fj(cj) affect the frequency- dependent occupancies Shj(uj) of all the receptors. 

Two receptors 

We first consider two identical receptors (Fig. 1) for which Eq. [16] simplifies to two coupled 
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linear equations 



— iuSni 



- n x 

n 2 ) 

- n 2 



2n 2 D 



- (fc+Ci + fc. 
5h 2 



A-nD |fi 



2it 2 D 



+ - 

- x 2 \ 

(k + c 2 + k. 



8h\ 



+ 



1 ~ Wi)Ci 

<5n 2 

knT 



8F X 



(17) 



5F 2 



47rD |fi — x 2 

In matrix form these equations can be written in terms of a complex susceptibility x as 



(18) 



5Fi 
<5F 2 



X 



5ni 



fc+(l— ni)ci 
—to; 
4nDrc2 



4irDrci 
fe+(l— n2)52 



<5n 2 



(19) 



where r = |xi — x 2 | and Ej 
to changes in the occupancies Sfii 



— . The inverse susceptibility x 1 relates changes in 5Fi 



2 ttDs 



i.e. 



Sfii 
8n 2 



5F, 



ft ffi w 



(20) 



According to the fluctuation-dissipation theorem [24|, the noise power spectra of the 



occupancies S, 



(5fij(uj)8hJ(uj)} can be calculated in the low-frequency limit from 



the above deterministic response functions, obtained by inverting the matrix x i n Eq. 



S nini (u —> 0) 



S I, 



ll-> 



uJ 



0) 



0) 



0) 



r 2k B T 
hm 1m 



r 1k B T 
hm 1m 



v 1k B T 
hm 1m 



v 2k B T 
hm 1m 

u>— >0 UJ 



dhA 2ife + (l-fii)ci(H-Ei) 
WJ (Mi + M 2 

= 2((5n 1 ) 2 ) tot r Cl (l + E 1 ) 

dn 2 \ 2fc+(l - n 2 )c 2 (l + Eg) 
dF 2 ) (k + c 2 + k.) 2 

= 2((5n 2 ) 2 ) tot r C2 (l + E 2 ) 
<9nA k\(l - ni)(l - n 2 )c 2 

ap 2 y 27iDr(k + c 1 + k_)(k + c 2 + k_) 

tot 

27rZ?rci 

<9ri 2 \ fc^l - ni)(l - n 2 )ci 

9^1/ 27r J Dr(fc + ci + kJ){k+c 2 + fc_) 
((^ 1 ) 2 ) tot ((5n 2 ) 2 ) tot 
27rDrc 2 



(21) 



(22) 



(23) 



(24) 



where ((<5n,,) 2 ) to t = «j(l — is the total variance in each occupancy and t Cj = (k+Cj + k-)^ 1 
is the correlation time for receptor occupancy in the absence of rebinding. The correlation 
time r c is modified when rebinding is included via coupling to particle diffusion so that the 
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effective correlation time for receptor occupancy increases to t Cj (1 + £j) as shown by the 
above equations. For an equilibrium system with time-reversal symmetry, the power spectra 
are necessarily symmetric in the zero-frequency limit, i.e. S nin2 (0) = S n2m (0). However, 
Snin 2 (0) (Eq. [23]) and S n2ni (0) (Eq. [21J are different depending on t\ and c 2 , respectively 
Hence, our approach based on the fluctuation-dissipation theorem is only valid for shallow 
gradients with c\ ~ c 2 . 

To determine the uncertainty of gradient sensing, we need to know the uncertainty in 
measuring a concentration difference. This uncertainty can be obtained from the noise 
power spectra for ligand binding and unbinding (Eqs. |2"TH!M]) . Starting with the expression 
n = c/(c + Kd) for the average occupancy in terms of the average concentration, where 
Kd = k^/k + , we obtain the following relation 



_ jc.+Kp) 2 _ K D _ 



which expresses how a fluctuation in the occupancy of receptor j affects the cell's best esti- 
mate of the local concentration. To obtain the uncertainty of the gradient measurement, we 
require the variance of the estimated concentration difference ([S(ci — C2)] 2 ), which we cal- 
culate using Eq. [25] and substituting the previously calculated noise power spectra S njni (0). 
For an averaging time r, we obtain for the variance of the inferred concentration difference 
between two receptors separated by a distance r: 

(Wd - c 2 )l) = ^ + KD)4 S nini (0) + (C ~ 2 +f D) \ 2W2 (0) (26) 



Kir 



2c, 2c 2 + £ L Wl_j_ ) (2?) 



k + (l — ni)r /c + (l — ri2)r itDt \s 2r 

Assuming that the true gradient is shallow, c\ ~ C2, to estimate the variance in Eq. [271 we 
set Ci = c 2 = c and n\ = n 2 = n for identical receptors in equilibrium, and obtain 

<[^-^)]?) = F7^V + ^ L ( 1 -f)' < 28 » 

k + {± — n)T itDt \s It) 

Two rings of receptors 



Following Bialek and Setayeshgar [24j . we next consider the case where a cell is equipped 
with two rings of m receptors, one at each pole. We further assume that the cell length is 
large enough so each receptor in Ring 1 is at nearly the same distance r from each receptor 
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in Ring 2 (Fig. 2). The equation governing receptor j in Ring 1 is 



Mi- 
+k + (i 

kjl 



2ttL)s 



+ 



4vrD 
ni)ci 



(/c+ci + k. 



Ring! 



Ringl \ Ring 2 ,^(2), x 

y, 6n]'{u) t ^ £n} >{uj) 

. Irt \ x 3 1=1 



(29) 
(30) 
(31) 



Summing over all receptors of Ring 1, and using the notation m8N^ = YlT=i Sn^ (u) and 



5F (i) = y™ =1 6FV(u), we obtain for Ring 1 



rnksT 



k + (l - nAci 



k + c\ + k- — iou[l + Si 



fe + (l - 
AttD 



$(Dl I SN W 



iujmkBT 
AixDr 



5N m = 5F (D 7 

(32) 

jj^x \x 1 -x \ '> w ^ n a similar expression for Ring 2. These equations can be 
solved for the uncertainty in measuring the concentration difference between the rings in a 
similar manner to the two-receptor case. For an averaging time r and $ = $W = & 2 \ one 
obtains for the uncertainty in measuring a gradient using two rings of m receptors each 



where = JJi 



([S(ci - c 2 )f T ) 



2ci 



+ 



2c, 



mk + (l — rii)r mk + (l — n 2 )r 



c\ + c 2 / 1 
vrnxDr \ s 



2" 



1 

2r 



(33) 



Assuming again that the true gradient is shallow, c\ ~ C2, to estimate the variance in Eq. 
531 we set Ci = C2 = Co and ni = — n for identical receptors in equilibrium, which results 
in 

2\ = 4c ° + 2c /I + £ _ J_ 
m/c_i_(l — n)r vhtxDt \s 2 2r 



(34) 



III. RESULTS AND DISCUSSION 

Many types of cells are known to measure spatial chemical gradients directly with high 
accuracy. In particular, Dictyostelium discoideum is known to measure extremely shallow 
gradients of cAMP important for fruiting body formation 0, @, 0] and Saccharomyces cere- 
visiae (budding yeast) detects shallow gradients of a mating pheromone {26|. Direct spatial 
sensing of gradients is also performed by cells of the immune system including neutrophils 
and lymphocytes jl], as well as by the large marine bacterium Thiovulum majus [9J. The 
question arises what are the limits of the accuracy of gradient sensing set by chemical dif- 
fusion? Recently, we derived fundamental physical limits for gradient sensing 



23 



using as 

model cells a perfectly absorbing sphere and a perfectly monitoring sphere (llj. We found 
that a perfectly absorbing sphere is superior to a perfectly monitoring sphere for both con- 
centration and gradient sensing since the perfectly absorbing sphere avoids the noise due 
to remeasuring previously detected particles (Table I). Consequently, our results for the 
perfectly absorbing sphere represent the true fundamental limits of both concentration and 
gradient sensing by cells. 

Our models of the absorbing and the monitoring spheres neglect all biochemical 
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reactions, such as particle-receptor binding and downstream signaling events, which 
might further increase measurement uncertainty. To study the effects of particle-receptor 
binding, we extended a formalism for the uncertainty of concentration sensing, recently 



developed by Bialek and Setayeshgar [24|, to gradient sensing. This formalism uses the 
fluctuation-dissipation theorem to infer the fluctuations of the receptor occupancy (and 
hence the accuracy of concentration sensing) from the linear response of the average 
receptor occupancy to changes in receptor binding free energies. The effect of particle 
rebinding is included by coupling particle-receptor binding to the diffusion equation (24j . 
leading to correlations in time among the receptors. 

Single receptor for concentration sensing 

It is instructive to first review the result for a single receptor without and with coupling 
to particle diffusion (corresponding, respectively, to preventing and allowing rebinding of 
already measured particles) [13]. The uncertainty of sensing concentration Co without re- 
binding (i.e. assuming that upon unbinding the particle is removed from the system) is 
given by 

where k + is the rate constant for binding, n is the average receptor occupancy, and r is the 
averaging time. The right hand side of Eq. [35] is obtained for diffusion-limited binding, in 
which case k + (l — n) — > AttDs with D the diffusion constant and s the receptor dimension. 
This is the well-known fundamental limit derived by Berg and Purcell jll|. In contrast, the 
uncertainty of concentration sensin g in cluding particle diffusion and possible rebinding is 



given by Bialek and Setayeshg [24 



/(x \2\ = 2 ( X + S ) c o = 2c o c o 3c 

U >t) k+(l-n)r h+(l - n)r itDst ^ 2itDst' 1 ' 

where £ = k + (l — n)/(2irDs). Comparison of Eq. [35] and Eq. [361 shows that the uncertainty 
of concentration sensing by a single receptor is larger by the term cq/(tiDst) when allowing 
for rebinding of already measured particles. For the minimum uncertainty case set by 
diffusion-limited binding, and given by the right hand side of Eq. [33 this additional term 
simply leads to a factor of 3 increase to the fundamental limit derived by Berg and Purcell 
(Eq. [35]). 

Two receptors for gradient sensing 

As derived in Methods (see Eq. |2~81) . we find that the uncertainty of gradient measurement 
is given by 

2 2 4c , 2cq fl 1 



K+(l — n)r z r nVr z T \s 2r J 

As expected, the larger the receptor-receptor separation r, the smaller the uncertainty in 
the gradient, because of the larger "lever arm" between receptors. Note that the result in 
Eq. [23 for the uncertainty in the gradient is independent of the magnitude of the actual 
gradient, including the case when there is no real gradient present. For comparison, the 
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uncertainty of mean concentration measurement is 24 



Analogous to the single receptor case, Eq. [36], the first term in Eqs. [37] and [38] arises from 
particle-receptor binding kinetics, whereas the second term is due to diffusion and includes 
the effects of possible rebinding of already measured particles. Due to the proximity of the 
receptors, separated by distance r, a particle can unbind one receptor and subsequently 
rebind the other receptor because of diffusion (see Fig. 1, trajectory 2). Hence, there is an 
additional noise component, which actually improves the accuracy of gradient measurement 
(term ~ — l/(2r)) due to cancellation with the noise due to rebinding to the same receptor, 
but degrades the accuracy of mean concentration measurement (term ~ +l/(2r)) since 
rebinding to the other receptor can only increase noise in the estimate of the mean 
concentration. This correlated noise was recently also investigated with Monte Carlo 



simulations 27 . 



Two rings of receptors for gradient sensing 

We next consider two rings of receptors, parallel to one another at opposite cell ends a 
distance r apart (Fig. 2). As derived in Methods (see Eq. [311) . we find that the uncertainty 
of gradient measurement is given by 

<«c: - *)]»/f = . ■ *> ■ , + (i + | - f) , (39) 

mk + {± — n)r z T m7rDr z T \s 2 2r J 

where m is the number of receptors per ring and $ is a geometric factor close to unity. For 



comparison, the uncertainty of mean concentration measurement is [24 



^ + ^ = ^tnfr + 2^57 (; + ! + h) ■ (40) 

The factor 1/m in Eqs. [39] and [40] reflects signal averaging by multiple receptors, which 
reduces the measurement uncertainty with respect to the case of two receptors. The 
possibility of rebinding to other receptors within the same ring leads to correlations among 
the signals, which are reflected in the extra term $/2 in the rebinding noise. 

Comparison with the perfect monitor and the perfect absorber models 

To make comparison to our results for the perfectly absorbing and monitoring spheres (Table 
I), which do not include particle-receptor kinetics, we replace k + (l — n)co by AttDscq for the 
minimum uncertainty case set by diffusion-limited binding. To specifically compare with the 
perfectly absorbing sphere, we neglect the second term in Eqs. 1391 and 1401 (thereby neglecting 
rebinding of particles) and obtain for gradient and concentration sensing 



(co/r) 2 TrDa'coT 
+ c 2 )/2] 2 ) = 1 
Co AkDcl'cqt 



(41) 
(42) 



11 



respectively. To specifically compare with the perfectly monitoring sphere we keep both 
terms in Eqs. EH] and HU] and obtain for gradient and concentration sensing 



(Wi-C2)]l)/r 2 

(c /r) 2 
(Wi + c 2 )/2]l) 



ixDcl'cqT 

1 

AirDa'coT 



3 + s 
3 + s I <3 



r 



(43) 
(44) 



The parameter a' = ms is the combined receptor dimension, ultimately limited by the cell 
dimension. Note that in Eqs. STJ and H3] for gradient sensing we normalized by (co/r) 2 , and 



in Eqs. 02] and HU for concentration sensing we normalized by Cq in order to use the same 
notation as Table I and Ref. (23| . 

As a result, for r > s, i.e. receptor separation larger than receptor size, the measurement 
uncertainty with rebinding (Eqs. [43] and 1441) is always larger than the measurement 
uncertainty without rebinding (Eqs. [4TJ and |42|) for both gradient and concentration 
sensing. Hence, the absorber is superior to the monitor even when receptor binding 
kinetics are explicitly included in line with our previous finding (Table I). Specifically, 
for diffusion-limited binding, the dominant effect of particle rebinding (Eqs. H31 and |4"4"|) 
is simply an increased numerical prefactor, also in line with our results for the perfect 
absorber and perfect monitor models. 

In conclusion, we found that the accuracy of concentration and gradient measurement 
without ligand rebinding is higher than the accuracy with rebinding, confirming the supe- 
riority of the absorber over the monitor (23| . Our model of two coaxial rings qualitatively 
resembles the polar clusters found abundantly in bacteria and archaea [28(. Hence, our 
model may be directly suitable for describing the concentration sensing by these organisms 
and possibly also for oxygen-gradient sensing by the bacterium Thiovulum majus [9|. Fur- 
thermore, a number of mechanistic models for gradient sensing and chemotaxis by eukaryotic 
cells have addressed the importan t q uestions of cell polarization, signal amplification, and 



adaptation [2 
gat ion [3 



3fl, 31, 32, 33, 34, 35 



cell movement of individual cells [3f 
as well as sensing of fluctuating concentrations [H, [2?], [39|, m. 



|37| . cell aggre- 
Our results on 



the accuracy of gradient sensing complement these models, and may ultimately help lead to 



a comprehensive description of eukaryotic chemotaxis [41 
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